Perturbations of the Yang-Mills field in the universe 



Wen ZhacEl 

Department of Physics, Zhejiang University of Technology, Hangzhou, 310014, China and 
School of Physics and Astronomy, Cardiff University, Cardiff, CF24 3AA, United Kingdom 

It has been suggested that the Yang-Mills (YM) field can be a kind of candidate for the inflationary 
field at high energy scales or the dark energy at very low energy scales, which can naturally give the 
equation of state — l<w<OorLi;<— 1. We discuss the zero order and first order Einstein equations 
and YM field kinetic energy equations of the free YM field models. From the zero order equations, 
we find that lo + 1 oc , from which it follows that the equation of state of YM field always 
goes to —1, independent of the initial conditions. By solving the first order Einstein equations and 
YM field equations, we find that in the YM field inflationary models, the scale-invariant primordial 
perturbation power spectrum cannot be generated. Therefore, only this kind of YM field is not 
enough to account for inflationary sources. However, as a kind of candidate of dark energy, the 
YM fleld has the 'sound speed' = — 1/3 < 0, which makes the perturbation (f> have a damping 
behavior at the large scale. This provides a way to distinguish the YM fleld dark energy models 
from other kinds of models. 
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I. INTRODUCTION 

It is well known that the physics of inflationary fields and 
dark energy are two of the most important problems in mod- 
eling cosmology. They all need a kind of matter with neg- 
ative pressure p ~ —p (Kolb & Turner [l^l; Bennett et al. 
[3|; Riess et al. Tegmark et al. [U; Feng et al. [T^l), 

where p and p are the pressure and energy density of matter 
respectively. People always describe them with a scalar field, 
which can naturally give an equation of state of — 1 < a; < 
(Kolb & Turner Wetterich [H; Bharat & Peebles 
In particular, when the potential of the scalar field is domi- 
nant, Lu will go to —1. Thus, the expansion of the Universe is 
close to the de Sitter expansion. Recently, the observations 
of the cosmic microwave background radiation (CMB) tem- 
perature and polarization anisotropies by the Wilkinson Mi- 
crowave Anisotropy Probe (WMAP) shows that the spectral 
index of the primordial density perturbation Us — 1.20lg jj 
at wavenumber k=0.002Mpc~^ (Bennett et al. Obser- 
vations of Type la Supernova (SNela), CMB and large scale 
structure (LSS) (Bennett et al. [J; Riess et al. [S^l; Tegmark 
et al. [21I I') also suggest that the equation of state of the dark 
energy may he lo < —1 (Corasaniti et al.Q]). These are all 
very difficult to obtain from single scalar field models. So it 
is necessary to look for a new candidate for the infiationary 
fleld and dark energy. 

Recently, a number of authors have considered using a vec- 
tor fleld as the candidate of the inflationary fleld or dark en- 
ergy (Ratra We have advised the effective YM conden- 
sate (Zhang [2I]; Zhang Zhao & Zhang [13]; Zhang et al. 
[2^) as a kind of candidate, which can be used to describe 
the inflation at high energy scales and dark energy at very 
low energy scales. In our models, a quantum effective YM 
condensate is used as the source of inflation or dark energy, 
instead of a scalar field. This model has the desired inter- 
esting feature: the YM field is an indispensable cornerstone 
to any particle physics model with interactions mediated by 
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gauge bosons, so it can be incorporated into a sensible unified 
theory of particle physics. Besides, the equation of state of 
this field is different from that of general matter as well as 
scalar fields, and the state uj < —1 can be naturally realized. 

In this paper, we shall discuss the evolution of the equation 
of state of the YM field and cosmic perturbations by solving 
the zero and first order Einstein equations and kinetic en- 
ergy equations. From the zero order Einstein equations, we 
find that the YM fleld can easily give a state of homogeneity 
and isotropy, and it can naturally give an equation of state 
LO < —1 01 LO > —1. We also find that t^; -I- 1 cx from 
the zero order equations. It follows that lo naturally goes to 
— 1 with the expansion of the Universe, independent of the 
initial condition. By considering the evolution of cosmic per- 
turbations, we investigate the first order Einstein equations 
and kinetic energy equations. In the simplest condition with 
only an electric field, we find that this YM field has the sound 
speed Cg = — 1/3 < 0, which is very different from scalar field 
models. We also find that a scale-invariant primordial per- 
turbation power spectrum cannot be generated, which shows 
that an electric YM field alone cannot be a candidate for the 
infiationary field. However, as a candidate of dark energy, 
YM fleld makes the cosmic fluctuation (p have a damping at 
large scales. This is helpful for answering large scale damping 
of the CMB anisotropy power spectrum. 



II. THE QUANTUM EFFECTIVE YANG-MILLS 
FIELD 

In the quantum effective YM field dominated Universe, the 
effective YM fleld Lagrangian is given by (Adler 01) 

C,s ^hF\n\F/eK'\, (1) 

where k is the renormalization scale with dimension of 
squared mass, F = —^F^^,F'^'"^ plays the role of the order 
parameter of the YM fleld. The Callan-Symanzik coefflcient 
b = llA''/247r^ for SU{N) when the fermion's contribution 
is neglected. For the gauge group SU{2) considered in this 
paper, one has 6 = 2 x ll/247r^. For the case of SU{3), the 
effective Lagrangian in Eq.([T| leads to a phenomenological 



2 



description of the asymptotic freedom for the quarks inside 
hadrons (Adler [11]). It should be noted that the SU{2) YM 
field is introduced here as a model for a cosmic infiationary 
field or dark energy, and it may not be directly identified as 
QCD gluon fields, nor the weak-electromagnetic unification 
gauge fields. 

An explanation can be given for the form in Eq.([T} as an ef- 
fective Lagrangian up to a 1-loop quantum correction (Pagels 
& Tomboulis 111]; Adler 01). A classical SU{N) YM field 
Lagrangian is £ = t^P^ where go is the bare coupling con- 
stant. As is known, when the 1-loop quantum corrections 
are included, the bare coupling go will be replaced by a run- 
ning g as used in the following (Gross&Wilczez Pagels 
& Tomboulis 0; Adler Q), gl ^ g" = = 
bin(fcVfc^) ' ^^^'^^ ^ '® momentum transfer and fco is the 
energy scale. To build up an effective theory (Pagels & 
Tomboulis To*]; Adler T]), one may just replace the momen- 
tum transfer fc^ by the field strength F in the following man- 
ner, ln(fc^/fco) — > 21n \FleK^\, yielding Eq.([lJ. The expression 
of is based on the renormalization group estimates. As 
emphasized in references (Adler Adler & Piran 01), this 
estimate is formally valid whenever the running coupling g^ 
is small in magnitude, which is true both when F j 2> 1 
(given 3^ is small and positive) and when F j ^ 1 (given 
g^ is small and negative). 

The attractive features of this efi'ective YM model include 
the gauge invariance, the Lorentz invariance, the correct trace 
anomaly, and the asymptotic freedom (Pagels & Tomboulis 
[r^V With the logarithmic dependence on field strength, /Icff 
has a form similar to the Coleman- Weinberg scalar effective 
potential (Coleman & Weinberg [^), and the Parker-Raval 
effective gravity Lagrangian (Parker & Raval [l3|)- The di- 
electric constant is defined by e = 2dCaS / OF , and in the 
1-loop order it is given by 



h\n\F/K 



(2) 



As analyzed in (Adler 0|), the 1-loop model is a universal, 
leading semi-classical approximation. Thus, depending on 
whether the field strength |F| > or \F\ < , the YM 
condensate belongs to the family of forms whose dielectric 
constant e can be positive or negative. The properties men- 
tioned above are still true even if 2-loop order corrections are 
taken into account, an essential feature of the effective model 
(Adler 1]; Adler & Piran [|). 

It is straightforward to extend the model to the expanding 
Robertson- Walker (R-W) spacetime. For simplicity we shall 
work in a spatially fiat R-W spacetime with a metric da^ = 
a^(r)(dr'^ — ^ijdx'dx-'), where we have set the speed of light 
c — 1, "/ij — 5j denoting background space is flat, and r = 
J{ao/a)dt is the conformal time. The dominant matter is 
assumed to be the quantum YM condensate, whose effective 
action is 5* = / Cejf «*(■'") d'^x, and the Lagrangian Ceff is 
defined in Eq.([l}. By variation of 5* with respect to g'^'^ , one 
obtains the energy-momentum tensor. 



Tm^ = -g^^-Fln\F/eK'\+€F;:,F^ 



(3) 



where the energy-momentum tensor is the sum of 3 energy- 
momentum tensors of vectors, T^^ — X^a^M"- order to 
keep the total energy-momentum tensor homogeneous and 
isotropic, we assume that the gauge fields are only functions 



of time t, and = |(TaA°(i) (Zhao&Zhang YM field 



tensors are defined as usual: 



C Ab AC 



(4) 



This tensor can be written in the form with electric and mag- 
netic fields as 
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From the definition in Eq.Q, we can find that Ef = E2 = Ej 
and Bl = Bl ^ Bl. Thus F has a simple form F = E^ 



where E' = E ti E? and B^ = Ylt 



B 



Inserting (|4]) in 



([3]), we can obtain the energy density and pressure of the YM 
field given by (Zhao&Zhang, [U 



P = 



\<E^ + B^) 
\<E^ + B') 



+ \h(E^ - B^) 



+ -^KB' 



■E' 



and 



p + p=-e(E' +B'). 



(6) 
(7) 

(8) 



Eq.Q follows a conclusion of this study: the weak energy 
condition can be violated: p -I- p < 0, by the effective YM 
condensate matter in a family of quantum states with the 
negative dielectric constant e < 0. The whole range of allowed 
values of F is divided into two domains with |_F| > and 
< K^, respectively. In the domain where \F\ > k^, one 
always has e > 0, so the WEC is still satisfied. The other 



domain is 
violated. 



F\ < k'' in which e < 0, so that the WEC IS now 



III. THE ZERO ORDER EQUATIONS 
A. The Zero Order Einstein Equations 



Let us firstly investigate the Friedmann equations, which 
can be written as 



8ttG 



(p + 3p). 



where "dot" denotes d/dt. If we consider the simplest case 
with only an "electric" field, as used in the previous works 
(Zhang et al. [2^; Zhao&Zhang [IJ), the Friedmann equa- 
tions can be reduced to 



(9) 



dp _ 4/3p 
da^ 3 ' 
where we have defined fH = e/b. If |/3| 1, one gets 

a^P = constant, (10) 

which follows the relation of the equation of state : 



Lu + 1 (X a 



(11) 



From this relation, we find that uu will run to the critical 
condition with uj — —1 with the expansion of the Universe. 
In the next subsectioin, we will find that from the zero order 
YM field kinetic equation, we can get the same result, which 
means that the YM field equation and the Einstein equations 
are self-rational. 
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B. The Zero Order Yang-Mills Field Kinetic 
Energy Equations 

By variation of S with respect to A^, one obtains the ef- 
fective YM equations (Zhao&Zhang [23]) 

9p(a*6 F""") - ,r'"'Al{a*e F""") = 0, (12) 

which can be simplified as (Zhao&Zhang, [13] ), 

a^Ee = constant. (13) 

If the Universe is dominated by an "electric" YM field, the 
kinetic energy equation follows as 

a^lSe'^^^ = constant, (14) 

When |/3| <C 1, then this equation becomes 

a^P{l + 13/2) ~ 0^13 = constant. 

which also follows as 

Lj + loca'^. (15) 

This is exactly the same as that of the zero order Einstein 
equations. 

IV. THE FIRST ORDER EQUATIONS 

A. The First Order Einstein Equations 

As a kind of candidate of inflationary fleld or dark energy, 
it is very important to study the evolution of perturbations 
in the YM field and cosmic fluctuations. In this section, let 
us consider the flat R-W metric with the scalar perturbation 
in the conformal Newtonian gauge 

ds^ = a^(r)[(l + 2(t>)dT^ - (1 - 2i!)'yijdx'dx^]. (16) 

The gauge- invariant metric perturbation tp is the Newtonian 
potential and (j) is the perturbation to the intrinsic spatial cur- 
vature. The first order Einstein equations become (Mukhanov 
et al. [Hi): 

-Sn{n(t> + ip') ^ 'iTTGa'^STll, (17) 

{n(p + ip'),, = atvGo^st^, (18) 

= -AivGa^STl (19) 

where Tl = a' /a, D = <?!>— V' and the prime denotes d/dr. From 
the definition of the energy-momentum tensor Q and the 
metric of (|16|l . one can get the first order energy- momentum 
tensor: 

STq = {e~b){B5B~E5E)+2eE5E+2(f){e~b)B^+ ^ ~^ ^ Se, 

^ (20) 

-5Ti = {e-h){-B5B + E5E)-2(t>{e~h)B'^ 

+ €[{2B2SB2 + 2Bs5B3 - 2E^SEi) + 4(B| + bI)<I>] 

+ liE^-B^)Se+^-^^6e, (21) 



and others are aU zero, where E5E — E1SE1 + E25E2 + E35E3, 
and similar for BSB. 

To obtain the gauge-invariant equations of motion for cos- 
mological perturbations in a Universe dominated by this kind 
of YM field, we insert the general equations of the energy- 
momentum tensor into (JT^. First of all, from the i — j [i ^ j) 
equation it follows that we can set tf) = tp, since 5Tj = 
{i 7^ j). Substituting the energy-momentum tensor STI^ into 
the general equations and setting tp — (j), we find: 

-3HiH(t> + 0')+V'0 

= 4TvGa'^{{£~b){BSB ESE) + 2eESE 

+ 2^ie~b)B' + ?-^5e}, (22) 

(H^-f 0').« = 0, (23) 

[{2H' + n^)ct) + H(t>' + (l>" + 2n<t)'\ 

= A-KGa{{-B5B + ESE){e - 6) - 2(t>B'^{e - h) 

+ y^4^B5B -2E5E + S.(t>B^] + ^ 5e}. (24) 

where Se = {2ESE - 2B5B - AB'^ (t>) / [E'^ - B^), the relation 
between 5E and 5B depends on the YM field equations as 
below. 



B. The First order kinetic equations of the YM 
field 

The metric as before with ijj = (f> is 

ds^ = a{T)[(l + 2(t>)dt^ - (1 - 2(l))j^jdx'dx^], (25) 

then ^y~g = a'^{l — 2(f>). The equation kinetic energy equation 
is 

d^{a\l - 2<^)e F""") - A^ia" (l - 2(f>)e F"'"') = 0, (26) 
from which, one can get the first order perturbation equations: 
d^,[a'^{edr5A + 5edrA-2e(f)drA)] = 0, (^i = 0,1,2,3) (27) 

^^ [5eA^ + 2e(f)A^ + 2eA5A] = 0, (i = 1, 2, 3) (28) 
which also can be written as 

di{BS€ + 2e(j>B + 2eSB) = 0, (29) 

^^{e5E + 5eE -2e(fiE) ^0, (30) 
dr [a^ {eSE + 5eE - 2ect>E)] = 0. (31) 

From these equations, we can immediately get the simple re- 
lation between SB and 5E: 

a' {eSE + SeE - 2e(l>E) = constant. (32) 

which is useful when solving the first order Einstein equations 
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C. The Solution of the Perturbations. 

Here we only discuss the simplest case with B = 0, the first 
order energy-momentum tensor becomes very simply 



STo = (e + 2b)E5E, 

- ST^ = (e - 2b)E5E/3, 
ST^ = -2e[Ei{5B3 - dB2)] = 



0, 



(33) 

(34) 
(35) 



where we have used Se = 2bESE/E^ when B = 0. From 



Ea. (l33|l - (|35|1 . we find that ST^f is independent of the metric 
perturbation of (j>, which is different from the scalar field mod- 
els (Weller & Lewis 22;]; Armendariz-Picon et al. DeDeo 
et al. ^). The first order Einstein equations become 



- 3H(H(t) + </>') 



2b), 



(36) 



[{2H' + 



{rL4> + 4>'),r^Q, (37) 

+ m<j>'] = A-KGaESEie - 26)/3. (38) 



From this we obtain the main equation, which describes the 
evolution of the metric perturbation (j> with time: 

(f)" + m[i - 7)(^' -I- + {2u' + n'^ - m'^j)(i) = o, (39) 

where 7 = ^^p^- From this equation, one can easily find that 
the evolution of (j> only depends on 7 and Ti., but not on the 
first order YM field kinetic equations. This is because we have 
only considerd the "electric" YM field. If we also consider the 
B components, the YM equation (38) will be used to relate 
dE and SB. 

For the inflationary field, the most important prediction 
is that the inflation can generate a scale-invariant primordial 
scalar perturbation power spectrum, which has been found in 
the CMB power spectrum and large scale structure. Now we 
shall firstly consider whether a Universe dominated by this 
YM field can also generate a scale-invariant spectrum using 
the general scalar infiationary models. 

If we consider the YM field as a kind of candidate of in- 
flationary field, the YM field should satisfy the following 
constraints: 1) Firstly, we should require that the infiation 
can exist, which means that the YM field has a state with 
uj < —1/3. Using the p and p in Eq.(12), one can easily get a 
constraint on the YM field: 



e < b. 



(40) 



2) The energy density of the YM field should be positive, 
which induces the second constraint on the YM field: 



e > 



(41) 



3) Due to equation ((39]), if we can define an adiabatic vacuum 
state at the very high frequency (fc —> 00), which requires that 
7 < 0, this would require the third constraint on the YM field: 



e > 26, or e < -26. 



(42) 



We find these three simple constraints on the YM field can- 
not be satisfied at the same time, so this model cannot gener- 
ate a scale-invariant primordial power spectrum as a general 
scalar field infiationary model. However, it is necessary to no- 
tice that this does not mean that the YM field cannot be the 
source of the infiation. Recently, some authors have discussed 
a kind of curvaton reheating mechanism in non-oscillatory in- 
fiationary models (Feng & Li ^). In this kind of model, the 



primordial spectrum and the reheating can be generated by 
the other curvaton field. So, although the YM field cannot 
generate a scale-invariant primordial spectrum, the YM field 
can also be the background field in the curvaton field infia- 
tionary models, which will be discussed in a future work. 

However, the YM field can be a good candidate of dark 
energy (Zhao & Zhang [13]; Zhang et al. [26]). Now let us 
discuss the evolution of the cosmic scalar fiuctuations in the 
YM field dark energy models. 

Neglecting anisotropic stress, the potential (f) evolves as 
Weller & Lewis ^ll; Ma & Bertschinger ^; Gorden & Hu 
M) 



+ 3H{1 



T-L^ + 2T-L' 



(43) 



where p — "^^Pi and p — "^^pi, which should include the 
contributions of the baryon, photon, neutron, cold dark mat- 
ter, and the dark energy. Here we consider the simplest case 
with only the YM field dark energy, and it has the equation of 
stateojde = —1. The 'sound speed' is defined by — 5p/5p. 
From the Eqs. (39) and (40), one finds that cl = -1/3. So 
then, the equation (|43p becomes: 



h" + 2n<j>' + - V^</> -I- 2H'(t) = 



which is the same as equation 
u = a(f), one gets 

" n 



(44) 

with 7 = 1/3. Defining 



(45) 



where Uk is the Fourier component with wavenumber k, we 
have used 00c 1/t, which has the solution of 



Uk oc . 



For the large wavelength case (fcr ^ 1), we have 
<f) oc 

and when fcr 3> 1, one has the growing solution 



u oc e , and 



a e n/3 . 



(46) 
(47) 

(48) 



We find in this universe, the evolution of the cosmic fiuctua- 
tion 4> is very different from that in the scalar field dark energy 
models. For the fiuctuation with large wavelength kr <^ 1, it 
is always damped with the expansion of the Universe. How- 
ever for the fiuctuation with small wavelength, it has a rapid 
growth with time, which is because of the negative sound 
speed Cs = — 1/3 < 0. This should have an important effect 
on the large scale integrated Sachs- Wolfe effect of CMB. We 
should mention that the result in this subsection is only qual- 
itative since we have not considered other effects, such as the 
baryons, photons, neutrons, dark matter and so on, which are 
also important for the evolution of the cosmic fiuctuation (f), 
especially at small scales. 



V. CONCLUSION AND DISCUSSION 

In this paper, we have investigated the effective YM field 
as a kind of candidate for an infiationary field or dark energy. 
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From the zero order Einstein equations and YM field kinetic 
energy equations, we find that this field can naturally give an 
equation of state — 1 < a; < and ui < —1. This is one of the 
most important feature of this model. Also, we get tj + 1 oc 
a~^, which suggests that lj will go to —1 with the expansion 
of the Universe. This makes the Universe dominated by this 
field have a de Sitter expansion. 

When considering the evolution of the perturbations, we 
solved the first order Einstein equations and YM field kinetic 
energy equations. We find that the model with only an elec- 
tric field cannot generate a scale-invariant primordial scalar 
perturbation power spectrum. This means that a model that 
only used the YM field is not a good candidate for describ- 
ing the inflationary field. However, as a kind of candidate 
of dark energy, we obtained the equation of the YM pertur- 
bation, and found that this field is very different from the 



general scalar field dark energy models. This YM field has 
a negative sound speed, which makes the cosmic perturba- 
tions dampen at large scales. This is a very important source 
of the integrated Sachs- Wolfe effect of the CMB power spec- 
trum. This damping of </> at large scale may be helpful for 
answering the very small quadrupole problem of the CMB 
temperature anisotropy power spectrum. 

We should mention that, in this paper we have not con- 
sidered the possible interaction between the YM field dark 
energy and the other components, especially the dark matter 
(Zhang et al. ^^). We leave this topic to future work. 
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